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Abstract 

Dynamics and stability of average current-mode control of buck converters are analyzed by sampled- 
data and harmonic balance analyses. An exact sampled-data model is derived. A new continuous- 
time model "lifted" from the sampled-data model is also derived, and has frequency response 
matched with experimental data reported previously. Orbital stability is studied and it is found 
unrelated to the ripple size of the current-loop compensator output. An unstable window of the 
current-loop compensator pole is found by simulations, and it can be accurately predicted by 
sampled-data and harmonic balance analyses. A new S plot accurately predicting the subharmonic 
oscillation is proposed. The S plot assists pole assignment and shows the required ramp slope to 
avoid instability. 
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1 Introduction 



Subharmonic oscillation is a fast-scale instability commonly seen in DC-DC converters. The insta- 
bility is generally associated with period-doubling bifurcation (PDB) [1]. Sampled-data models [2HZ] 
can accurately predict PDB, which occurs when a sampled-data pole is at -1 in the complex plane. 
A sampled-data pole p, not negative real, can be mapped to a corresponding pole \n(p)/T in the av- 
erage model. However, there is no single corresponding pole which can be mapped from a negative 
real sampled-data pole [8]. That explains why some average models cannot predict PDB. 

In the past research [91-11 7j. frequency responses of the average models show good matches 
with experimental data. However, these experiments are made under specific conditions when the 
sampled-data poles are not negative real. The deficiency or inaccuracy of the average model will be 
pronounced if there exists a negative real sampled-data pole. This deficiency can be eliminated by 
mapping ("lifting") the negative real pole to a pair of complex poles in the right half plane (RHP) 
and a zero [8], which increases the system dimension. With an additional pole, the lifted model 
has additional degree of freedom to have matched frequency response. This approach is adopted in 
this paper. An exact sampled-data model is first derived, then lifted to a continuous-time model. 
The lifted model is more accurate than the average model, and can predict PDB. 

For a converter under peak current-mode control (PCMC), the current loop is one-dimensional. 
By considering the sampling effect and increasing the current loop dimension to two, improved av- 
erage models [T8l[T9] can accurately predict PDB. In essence, the accuracy is improved by "lifting" 
the sampled-data model to increase the system dimension. In contrast, any improvement in the 
average model without increasing the system dimension is deficient in predicting PDB. The oscilla- 
tion predicted by [20] is actually a Neimark bifurcation (slow-scale instability), not a PDB, because 
the predicted oscillation frequency is not half the switching frequency. A single continuous-time 
pole cannot cause system oscillation, while a single negative real sampled-data pole can. A pair of 
RHP continuous-time poles are required to cause oscillation. Without increased system dimension 
as in [18[ll9j. the average model generally cannot predict PDB. 

Previous works on average current-mode control (ACMC) are generally based on average models 
[91117] . which can be divided into two groups. The first group [91415] does not consider the sampling 
effect and is equivalent to state-space average models. The second group [THE] considers the 
sampling effect. Adding the sampling effect is effective for PCMC, but directly applying it to 
ACMC has been questioned |10j . Two reasons can be given. First, ACMC resembles voltage-mode 
control rather than PCMC [10] , and considering the sampling effect may be unnecessary. Second, 
The current loop in ACMC has more than one pole, not just one as in PCMC. Some poles may 
migrate and regroup with each other |12| I14| ] T5]. One cannot predict which pole is real negative to 
be lifted. It is not accurate as in [16] to lift just one particular pole. 

Harmonic balance modeling |21| is another approach to analyze switching converters and can 
predict PDB. The contributions of the paper are as follows. First, PDB does occur in ACMC. It even 
occurs in the circuit example in [IB], and such a PDB instability was reported around a decade ago 
in [22]. The occurrence of PDB can be predicted by the sampled-data and harmonic balance models. 
Second, it is shown that the average model has frequency responses matched with experimental 
data only in those conditions when the converter does not have negative real sampled-data poles. 
The average model without increased system dimension fails to predict PDB. Third, a lifted model 
is proposed. It has frequency responses matched with experimental data, and it accurately predicts 
PDB even if there is a negative real sampled-data pole. Fourth, a new stability criterion based on 
the harmonic balance model is derived. Fifth, a new S plot is proposed. It accurately predicts PDB 
and greatly assists converter designers to adjust some particular parameters to avoid PDB. It also 
assists pole assignment. 
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Figure 1: System diagram of an ACMC DC-DC converter. 

2 Operation and Instability of ACMC 

An ACMC DC-DC converter is shown in Fig. [TJ The source voltage is v s , and the output voltage 
is v . Let the switching frequency be f s = 1/T and let uj s = 2irf s . The inductor current %l is 
sensed by a resistor R s and compared with a control reference v c . The difference is amplified by a 
compensator, 

HM ~ 1W) = C + W + dfc) ( ' 

where if c is the compensator gain at low frequency. The mid- frequency gain is K c /oj z . The poles 
are 5 ~ and uj p . The compensator output y(t) is compared with a T-periodic ramp signal h(t), 
which has h(0) = and h(T~) = Vh (ramp amplitude). The nominal solution of a DC-DC converter 
is a T-periodic orbit, not an equilibrium point as depicted in the average model. Denote the orbit 
by x°(t), and let y°(t) = Cx°(t) + Du. A periodic orbit (like x°(t) or y°(t)) is orbitally stable if a 
state trajectory stays in the orbit after being perturbed. If the state trajectory eventually leaves 
the orbit, the orbit is unstable. Previous work [10] reported that the ripple size or slope of y affects 
stability. However, the following simulations show otherwise. Furthermore, an unstable window of 
uj p is found. 

Example 1. (Unstable window of oj p .) Consider an ACMC buck converter [16^ p. 114]. The 
system parameters are v s = 14 V, v Q = 5 V, v c = 0.5, f s = 50 kHz, L = 37.5 /iH, C = 380 fiF 
with equivalent series resistance (ESR) R c = 0.02 O, R = 1 Q, R s = 0.1 O, 14=1, h(d) = 50000, 
K c = 75506, ujp = 0A92u s , and uj z = 5652.9 rad/s. All of the parameters are the same as in [16] . 
The converter is actually unstable [22]. The sampled-data poles are -1.123, -0.045, 0.882, and 
0.9537. The pole -1.123 indicates that PDB occurs and the orbit x°(t) is unstable. When x°(t) is 
perturbed, the state trajectory will leave the periodic orbit and become erratic as shown in Fig. [2] 
and does not converge to a T-periodic orbit. Time simulation shows that the converter is unstable 
for 0.13 < ojp/uis < 0.56. This unstable window of uj p can be predicted by the sampled-data and 
harmonic balance analyses shown later. □ 

Example 2. (Stability is unrelated to ripple size.) In Example [H change L to 46.1 /iH, which 
does not change the dynamics substantially. The change of L is used to illustrate a slightly smaller 
unstable window of uj p . Compensator pole u> p is varied from 0.1w s to 0.8w s . The pole u> p adjusts 
how much high-order harmonics of the compensator output y is attenuated. The ripple size and 
slope of y increase as uj v increases. An unstable window of u p was found and reported in [22]. In the 
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Figure 3: Unstable T-periodic y (t), u p = 0.49^. 



window 0.18w,g < u p < 0.49w s , x°(t) and y°(t) are unstable. Outside the window, x°(t) and y°(t) 
are stable. Therefore, the stability is unrelated to the ripple size: for u p < 0.18o; s , the converter is 
stable with a small ripple size of x°(t) or y°(t); for 0.18w s < uj p < 0.49w s , the converter is unstable 
with a little larger ripple size; for u p > 0.49o; s , the converter is stable again with a much larger 
ripple size. 

Let u p = 0.49w s , for example. Subharmonic instability occurs, and the unstable T-periodic orbit 
(Fig. [3]) and the stable 2T-periodic orbit (Fig. U]) coexists. When the unstable T-periodic orbit is 
perturbed, it will lead to the stable 2T-periodic orbit. The two orbits have the same average duty 
cycle 0.357 and share the same average orbits. Based on the average model, they should have the 
same stability. However, their stabilities are different. For uj p > 0.49a; s , x°(t) and y°(t) are stable. 
Let Up = 0.8u s , for example, y°(t) is stable with a large ripple size and slope (Fig. [5]). Comparing 
Figs. [3] and one can see that the stability is unrelated to the ripple size. It is interesting to note 
that as Up varies from 0.1u s to 0.8w s , the average duty cycle remains the same, but the stability 
can be quite different. □ 
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Figure 4: Stable 2T-periodic orbit, cj p = 0.49u;, 
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Figure 5: Stable T-periodic y°(t) with Zarge ripple, u p = 0.8lu s . 



3 Exact General Sampled-Data Analysis 
3.1 Exact Switching Model 

The switching model [6] in Fig. [6] describes exactly the operation of a DC-DC converter. In the 
model, A lt A 2 G R JVxAr , £i,5 2 G R Arx2 , C,E U E 2 G R lxJV , and I? G R lx2 are constant matrices, 
where iV is the system dimension. For example, JV = 4 for a buck converter with a two-pole 
compensator. Within a clock period, the dynamics is switched between S\ and S 2 - Switching 
occurs when the ramp signal hit) intersects with the compensator output y := Cx + Du G R. The 
system is in Si when y[t) > h(t), and switches to S 2 at instants when y(t) < h(t). 

In the n-th cycle, let x n = x{nT). For high switching frequency, u = (v s ,v c )' G R 2xl can 
be considered constant within the cycle and denoted as u n = (v sn ,v rn )'. The short notation v sn , 
instead of v St7l , is used for brevity. This applies to other variables. 

Let the steady-state duty cycle be D and let d = DT. Confusion of notations for capacitance 
C and duty cycle D with the matrices C and D can be avoided from the context. Let x°(d~) = 
Aix°(d) + B\u and x°(d + ) = A2X°(d) + B 2 u denote the time derivative of x°(t) at t = d~ and d + , 
respectively. 
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Figure 6: Exact switching model for DC-DC converter. 



The periodic orbit x°(t) in the system of Fig. [6] corresponds to a fixed point x°(0) in the sampled- 
data dynamics. Using a hat A to denote small perturbations (e.g., x n = x n — x (0)). From [6], the 
linearized sampled-data dynamics is 

X n+ i = + Tu n = $X n + TiV sn + T 2 Vcn (2) 

where $ G R^*^ and T = [ri,r 2 ] G R^^ are 

V Cx°(d-) - h{d) ' V ; 

r = e^-^/^^^-^^W. 

Vo Cx°(d-)-h(d) 

(C [ d e Aia daB 1 + D))+ f * e A ^daB 2 (4) 
Jo Jo 

The periodic solution x°(t) is asymptotically orbitally stable if all of the eigenvalues of $ are 
inside the unit circle of the complex plane. PDB occurs when an eigenvalue is at -1. Saddle-node 
bifurcation (SNB) [T| occurs when an eigenvalue is at 1. 

Let E := (E\ + E 2 )/2, because the output voltage may be discontinuous. From ([2]), the control- 
to-output transfer function is 

T oc (z) = V -^\=E{zI-<S>)- l T 2 (5) 

v c {z) 

Given a transfer function in the z domain, T(z), its frequency response is T(e^ T ), valid up to half 
the switching frequency. 

3.2 Lifted Model with Improved Frequency Response 

In the sampled-data model, the switching frequency is high enough that v c is assumed constant 
within the cycle, which is the same as assuming a sampler with a zero-order-hold being placed after 
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v c . The sampled-data dynamics (J2j> can be converted back ("lifted") to a continuous-time dynamics 
x = Ax + Bu, with the following relationship [H[5l[8], 



$ r 
o 1 



exp( 



A B 




T) 



This conversion can be done by the d2c command in Matlab. A negative real pole, if any, in the 
sampled-data model is mapped to a pair of complex RHP poles and a zero, and the lifted model 
has increased system dimension [8]. The "lifted" control-to-output transfer function is 



E(sl - A)~ 1 B. 



(6) 



The "lifted" model has improved high frequency phase response because the effect of the assumption 
that v c being constant within the cycle becomes significant at high frequency. It will be shown (in 
Examples 6 and 7) that the "lifted" control-to-output transfer function has phase response matched 
with experimental data. The lifted model also predicts PDB. 

3.3 "S Plot" Accurately Predicts Pole Location and PDB 

Based on [H p. 46], suppose A is not an eigenvalue of e A 2{T~ d ) e Md^ then A is an eigenvalue of <J> if 
and only if the following boundary condition holds: 



y°(<T) + C(Xe- A ^ T - d ^e- Ald - J)" 1 ^ ^ - ) - x°(d + )) = h{d) 
The proof is as follows. Suppose A is not an eigenvalue of e A2 ^ T ~ d ^ e Ald , then 
det[A/ - $] = det[A/ - e MT-d) e A ld] det[/ + (A/ _ e A 2 (T-d) e A ld yi e A 2 (T-d) 



(7) 



x°(d- 



b°(d + ) 



Ce 



Aid] 



Cx (d~) - h(d) 

= det[A7 - e A ^ T ~ d h Md ][l + Ce Md {XI - e MT-d) e A ld) -i e A 2 {T-d) *°(d~) ~ i°(d+) 
1 JL 1 } y°(d~)-h(d) J 

det[A7 — <£] = requires that the last term (inside the second square brackets) of the last equation 
equals to zero, which leads to ©. 

Since the condition ([7]) is in terms of signal slopes (y°(d~), x°(d~), x°(d + ), and h(d)), the left 
side of (|7|) is called an "S plot" , S(X, D, oj p ), as a function of A, D = d/T, and ui p , for example. Then, 

PDB occurs when S(— 1, D,uj p ) = h(d). The converter is unstable with PDB if S(— 1, D,uj p ) > h(d). 
SNB occurs when S(l, D,u p ) = h(d). Although this paper focuses on the buck converter, the 
condition (|7|) is also applicable to other converters. 

Like the popular Bode plot, the S plot facilitates the converter design. For example, the 
intersection of the S plot (as a function of pole A) and h(d) shows the real pole location. The S plot 
can be used for pole assignment. A limitation of the pole assignment is that only real poles can be 
assigned. Choosing a particular ramp slope h(d) can assign the pole to some locations. Also, if the 
ramp slope h(d) is large enough such that h{d) > S(—1,D,uj p ), PDB is avoided. 
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4 Modeling of ACMC Buck Converter 
4.1 Sampled-Data Analysis 

For the buck converter, let x = (i^, vq, v e i,v e 2)', where is the inductor current, vq is the capacitor 
voltage, and v e \ and v e 2 are the states of the current-loop compensator. Then, 
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where p = R/(R + i? c ). A small nonzero 5 is chosen to make both A\ and A 2 invertible. Let 
B\ := , -B12] 5 B2 '■= [-B2IJ-B22] t° expand the matrices into two columns. The buck converter 
generally has A\ = A2, B 2 i = OnxX, and B12 = £?22- The boundary condition ((7|) becomes 

5(A, D, u; p ) := C[(I - e MT )-\e A ^ DT -/)+/ + (Ae" AlT - iy l ]B llVs = h{d) (8) 

If v is fixed, i> s = v a /D is used in the S plot. 

4.2 Average Analysis 

The linearized dynamics of the state-space average model for the buck converter is derived: 



x (A . v s B n C ,v s B n v A 
x = (Ai H — )cc + (— h S 12 )u c 



The control-to-output transfer function is 



T - {s) ~us)- E{sI ~ Al — vr ] { — +Bl2) 



0) 



(10) 



4.3 Harmonic Balance Analysis and the S Plots in Frequency Forms 

Let the current loop transfer function be T(s). Then, one has T(s) = H c (s)Gid(s) /Vh, where 

v s ((R + R c )Cs + l) 



G id (s) 



(R + R c )LCs 2 + (L + RR c C)s + R 



(11) 



is the duty-cycle-to-inductor-current transfer function. Define a transfer function G{s) = T(s)Vf l /v s 
From [2TJ, PDB occurs when 



2 Us / a Re[£((l - e^ kD )G(jco s k) - G(jco s (k - -)))] = h(d) 
fe=i 1 



(12) 
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One can prove that (|12p and ([8j) (with A = — 1) are equivalent but expressed in different forms, and 
the left side of (fT2|) is the S plot, S(— 1, D, ui p ). 

Generally G(juj s ) is low-pass, and using just one term in the summation in (|12|) . 

S(-l,D,u p ) « 2v s / fl Re[(l - ^ D )G(ju s ) - G{^-)) (13) 

Using © and (HH), for oo s » l/v 7 ^ and 1/RC7, 

cv 1 n \~/ 4w p^ ,K c R s v 

s( - 1 ' A ^ ) " ( 4^f ) ^l^ (14) 

For 5(— l,D,uj p ) > h(d), instability occurs with PDB. A higher value of S(— l,D,u p ) is likely to 
lead to PDB. Generally, the S plot as a function of uj p is A-shaped (see Fig. [8] for example), and 
has a maximum value of K c R s v d /uj z LttD at uj p = u s /2. To avoid PDB, the ramp slope h(d) should 
be higher than this maximum value. Therefore, a conservative (conservativeness because it is valid 
for all uj p ) design guideline to avoid PDB is 

K c R s v s K c R s v Q ■ 

-^JT = -^LD < nk{d) (15) 

which shows how different parameters affect stability and is useful for converter design. If the 
condition (|15p is not met, there may exist an unstable window of uj p . Although the unstable window 
of u! p predicted by (fPo"|) may be a little different from the exact condition (Tl2"j) . the approximate 
condition (I15f) at least predicts a characteristic of an unstable window of uj p . No existing stability 
criterion has been reported to predict the unstable window of d p . 

It is reported in [15 J that an increase of the PWM modulator gain 1/Vh, the mid- frequency 
current loop gain K c /u z , or v s will cause a shift of resonance frequency toward / s /2. It is also 
reported in [15] that the load resistance R does not cause such a shift. These observations agree 
with (fT5j) . 

Note that (|12|) is an exact condition, (|13|) and (|14p are approximate, and (|15p is conservative. 
Also note that (|14p is w p -dependent, and such a dependence is removed in (|15p . Using which one of 
the conditions f|12|) - (|15[) depends on how much accuracy is required. A design guideline proposed 
in [12] suggests that the value of oj p is chosen between 0.33o; s and 0.5w s . However, this range is 
close to the unstable window of u p . 

Another stability guideline proposed in [10] suggests 



K c R s v s . 2 1 ■ 

— < mm -, —\h{d) (16) 

oj z L 1 — D D 

For a small D, the term 2/(1 — D) dominates and has a value between 2 and 3. For a large D, 
the term 1/D dominates and has a value between 1 and 3. Either way, the right side of (|16p has a 
value of at most 3. Therefore, (|16p is more conservative than (|15p . The stability criterion can be 
relaxed as in (fT5]) to achieve better dynamics. 



5 Illustrative Examples 

5.1 Many Useful Usages of the S Plot 

Example 3. [The S plot predicts the pole location and shows the required ramp slope to avoid 
PDB.] Consider again Example 1. In Fig. \7\ the intersection of the S plot (as a function of A) and 



10 



x 10 4 

16 
14 
12 
10 

a. 
o 

tn a 

6 
4 
2 


-2 -1.5 -1 -0.5 

PoleJ. 

Figure 7: The intersection of S(X, 0.357, 0.492w s ) (solid line) and h(d) (dashed line) shows two 
poles at -1.123 and -0.045. 

h(d) shows the two negative real poles at -1.123 and -0.045. To avoid PDB, the S plot shows that 
the required ramp slope is 62000 (corresponding to Vh = 1-24), which shifts the two poles to -0.999 
and -0.051. To shift the poles from being negative real, the S plot also shows that the required 
ramp slope is 150000. With this slope, the sampled-data poles are now either complex or positive, 
and they are —0.224 ± 0.029j, 0.872, and 0.957, agreed with the prediction of the S plot. 

[The S plot predicts the unstable window.] In Fig. [81 the intersection of the S plot (as a function 
of oj p ) and h(d) shows exactly the unstable window of uj p : 0.13 < oj p /uj s < 0.56. In ACMC, the 
output voltage is regulated to be around v c , then D is varied if v s is varied. The S plot as a 
function of D shown in Fig. [9] shows that the converter is unstable for all D. In both Figs. 
and El a circle symbol indicates the circuit parameters (e.g., D = 0.357 and ui p /oj s = 0.492) used 
in Example [TJ The circle symbol is above the ramp slope line at 50000, indicating instability as 
reported in Example [TJ Both of the figures show that if the ramp slope increases to 62000 (or 
Vh = 1-24), the converter is stabilized. 

With h(d) = 62000, the sampled-data poles are -0.999, -0.051, 0.881, and 0.9537, indicating 
stability as predicted. In Fig. [9j draw a horizontal line at h(d) = 62000. The line intersects with 
the S plot at D = 0.35 and 0.72. Therefore, this larger ramp slope leads to a larger stable operating 
range: 0.35 < D < 0.72. The converter in Example [TJ has D = 0.357 which is now within the stable 
range. 

In Fig. [HJ draw a horizontal line at h(d) = 62000. The line intersects with the S plot at 
ujp/uj s = 0.18 and 0.49. Therefore, this larger ramp slope also leads to a smaller unstable window: 
0.18 < ojp/u s < 0.49. The converter in Example [TJ has u p = 0.492cjs which is now outside the 
unstable range. This window is the same as in Example 2 (with increased L). This can be 
explained by f)15|) because increasing L or h(d) has the same effect to stabilize the converter. □ 

In summary, the S plot has many usages: 

1 . The S plot as a function of A (Fig. [71 for example) predicts the pole location and shows the 
required ramp slope to assign the real poles to some locations. 

2. The S plot as a function of uj p (Fig. [SJ for example) predicts the stable operating range of uj p 
and shows the required ramp slope to operate in a desired range of lo p without occurrence of 
PDB. 
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Figure 8: The intersection of <S(— 1, 0.357, lj p ) (solid line) and h(d) (dashed line) shows the unstable 
window of uj p G (0.13, 0.57)a; s . The circle, for the set of converter parameters in Example[H is above 
h(d) = 50000, implying instability in Example [TJ 
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Figure 9: S(— 1, D, 0.492o; s ) (solid line) is above h(d) (dashed line), showing instability for all D. 
The circle, for the set of converter parameters in Example [U is above h(d) = 50000, implying 
instability in Example [U 

3. The S plot as a function of D (Fig. [9l for example) predicts the stable operating range of D 
and shows the required ramp slope to operate in a desired rang of D without occurrence of 
PDB. 

If the S plot is represented as a function of another variable (K c , uj z , or L, for example), it can 
predict the stable operating range of such a variable and show the required ramp slope to avoid 
PDB. 



5.2 The Sample-Data Model Accurately Predicts the Unstable Window of u p 
Example 4. Consider again Example 1. The sampled-data model © is 

0.87528(z + 0.4034)(z - 0.8987)(z - 0.0255) 



(z + 1.123)(z - 0.9537)(z - 0.882)(z + 0.04509) 



(17) 
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Figure 10: Samp led- data pole trajectories for u p /oj s G (0.1,0.8). 

There are two negative real poles. One negative pole is less than -1 and is unstable. The lifted 
model ([!]) is 

29505Q + 5338)(s 2 + 302500s + 3.812 x 10 10 )(s 2 + 124700s + 2.981 x 10 1Q ) 
(s + 6276)(s + 2372)(s 2 + 309900s + 4.868 x 10 10 )(s 2 - 11620s + 2.471 x 10 10 ) ^ ' 

The unstable negative sampled-data poles are lifted to become a pair of unstable complex poles, 
accurately predicting the instability. The stable negative sampled-data poles are lifted to become 
a pair of stable complex poles. The lifted model now has six poles. 

Next, the effect of uj p is considered. The sampled-data pole trajectories for uj p /oj s 6 (0.1,0.8) 

— T 

are shown in Fig. [10J Two poles are fixed around 0.882, and 0.9537 (~ crc). A pole leaves the unit 
circle through -1 when u p = 0.13a; s , and enters the unit circle when u p = 0.56o; s . This explains 
exactly the unstable window of uj p discussed in Example [TJ Similar techniques to increase the 
system dimension to predict PDB is also adopted in |23,[24]. □ 

Example 5. {Average models do not predict PDB and the unstable window of u p .) In Example 
1, the average model (fT0|) is 

7320(s + 131600)(s + 2216000)(s + 5272) 
(s + 5945) (s + 2477) (s 2 + 149400s + 7.641 x 10 10 ) ^ ' 

with poles in the left half plane (LHP) and fails to predict the instability. 

By considering the sampling effect but without increasing the system dimension, another average 
model of [16] is 

3.6 x 10 14 (s + 131600) 
(s + 1962000)(s + 154700)(s + 12570)(s + 1243) ^ ' 

with poles still in LHP and fails to predict the instability. One explanation that these average 
models fail to predict PDB is that these average models do not lift the unstable sampled-data poles 
to increase the system dimension. 

Next, the effect of oj p is considered. The pole trajectories for u p /oj s G (0.1,0.8) of the average 
model (jlOp are are shown in Fig. [TTJ Two poles are fixed around -5945, and -2477, which can 
be mapped from the corresponding sampled-data poles 0.882 and 0.9537. For example, -5945~ 
ln(0.882)/T. All four poles are in LHP. This contradicts the unstable window in Example 1. □ 



13 



-1 - 

-2 - 




-14 -12 -10 -8 -6 -4 -2 

X 10* 

Figure 11: Average model pole trajectories for lo p /uj s £ (0.1,0.8). 

5.3 The Lifted Model Predicts PDB and Has Matched Frequency Response 

Example 6. Consider another ACMC buck converter from flUI p. 982]. The system parameters 
are (j s = 5V,u = 2V,« c = 0.279, f s = 180 kHz, L = 13 /iH, C = 750 fjF, R c = 5 m£l, R = 0.43 
U, R s = 0.06 n, 14=2.7, h(d) = 486000, K c = 98000, oj p = oj s , and oj z = 6723 rad/s. 
The sampled-data model (jSJ is 

0.078161(z - 0.2163)(z - 0.9654)(z + 0.08168) 

(21) 



(z - 0.003783)(z - 0.5155)(z - 0.9525)(z - 0.9861) 

All four poles are positive real. Since there are no negative real poles, the average model is expected 
to be accurate. The lifted model (jHJ) is 

8100(s + 1427000) (s + 259500) (s + 6332) 

(22) 



(s + 1004000) (s + 119300) (s + 8755) (s + 2528) 
There are still four poles, directly mapped from the sampled-data poles. The average model (|1Q[) is 

704(s + 17610000) (s + 266700) (s + 6294) 



(s + 989500) (s + 133800) (s + 8523) (s + 2573) 



(23) 



which has similar corresponding pole locations. The control-to-output frequency responses of the 
lifted and the average models are shown in Fig. [T2l both matching with the experimental data 
shown in Fig. 6 of [TO] . 

Here, the average model shows matched frequency responses. The good match in one condition 
does not necessarily mean that the good match occurs in all other conditions. The good match is 
due to the fact that here the sampled-data poles are not negative real. If one converter parameter 
(v s , for example, as shown next) is adjusted such that the converter has a negative real sampled-data 
pole, then the average model fails to predict the converter dynamics. 

The S plot (as a function of D) in Fig. [131 shows that PDB occurs if D < 0.065 (corresponding 
to v s = v /D = 30.84). The approximate S plot (|13p fits closely with the exact plot based on ([8]) or 
(|12p . The other approximate S plot (|14p. although not fitting with the exact plot, still accurately 
predicts the instability around D = 0.065. 
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Next, let v s = 30.84. The unstable time waveform is similar to Figs. 2 and 4, omitted to save 
space. The sampled-data model ([5]) is 

0.15520 + 1.692)0 - 0.9657)0 - 0.3773) 

(24) 



+ 1.0002)0 + 0.001935)0 - 0.9623)0 - 0.9835) 

which shows PDB with two negative real poles. One negative pole is unstable. Lifting these two 
poles to totally have six poles and five zeros, with the other two stable poles being not substantially 
changed, the lifted model ([6]) is 

137269Q + 1335000)0 + 258800)Q + 6284)(s 2 - 125700s + 2.522 x 10 11 ) 
(s 2 - 68s + 3.198 x 10 n )(s 2 + 2249000s + 1.584 x 10 12 )(s + 6920)0 + 2994) ^ ' 

which also has a pair of unstable poles and accurately predicts the instability. The average model 
CED is 

43430 + 266700) (s + 6294) (s + 17610000) 



2 + 1124000s + 8.62 x 10 n )(s + 7013)0 + 2963) 



(26) 



with poles in LHP and fails to predict the instability. The control-to-output frequency responses of 
the lifted and the average models are shown in Fig. [TH The lifted model has a peaking at half the 
switching frequency and shows the instability. The average model does not show such a peaking and 
fails to predict the instability. Therefore, even though the average model shows matched frequency 
response in one condition (v s = 5), it does not necessarily show accurate frequency response in 
another condition (v s = 30.84). In contrast, in either condition, the lifted model accurately predicts 
the dynamics. 

With v s = 30.84, the S plot as a function of lo p in Fig. [15] shows an unstable window for 
0.2 < ujp/uj s < 1. In this example, uj p equals to uj s and is inside this unstable window. To avoid the 
unstable window, one has two options: increase the ramp slope or operate with ui p being outside 
the window. 

First, consider the option of increasing the ramp slope, the S plot shows that the required slope 
is h(d) = 700000 (corresponding to Vh = 3.89). Let h(d) = 700000, the sampled-data poles are 
-0.652, -0.003, 0.962, and 0.984, indicating that PDB is avoided as expected. 

Next, consider the option of operating with uo p being outside the window. The S plot also shows 
that if uj p /uj s > 1 or uj p /uj s < 0.2, PDB is avoided by operating outside the unstable window. Let 
oj p = 0.15w s , for example, the sampled-data poles are — 0.457 ± 0.442j, 0.962, and 0.984, indicating 
that PDB is avoided as expected. □ 



Example 7. Consider another ACMC buck converter from [12]. The system parameters are 
v s = 3V,v = 2.25 V, v c = 0.5625, f s = 100 kHz, L = 20 fiR, C = 330 /jF, R c = 25 mO, R = 2 Q, 
R s = 0.5 n, V h =1.8, h(d) = 180000, K c = 11455, u p = 314940 » u,/2, and uj z = 11905 rad/s. 
The sampled-data model ([5]) is 

0.0474670 - 0.0746)0 - 0.9432)0 + 0.003245) 

2 - 1.564* + 0.6236)0 - 0.06766)0 ~ 0.9889) ^ ' 

All four poles are stable. Different from Example (6[ a pair of complex poles exists. This example 
is used to show that no matter how the poles migrate and regroup, both the sampled-data and the 
lifted models still accurately predict the dynamics. The lifted model (JSJ) is 

2739(s + 273900) (s + 209300) (s + 5847) 
(s 2 + 47230s + 758500000)(s + 269300)(s + 1119) ^ ' 
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Figure 12: Control-to-output frequency responses; solid line for the lifted model and dashed line 
for the average model, v s = 5. 
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Figure 13: The intersection of S{— 1,D,lj s ) (solid line) and h(d) (dashed line) shows instability for 
D < 0.065; dash-dotted line for (fT3|) and dotted line for (fT4"]) are approximate S plots. 



There are still four poles (with complex poles), mapped from the sampled-data poles. The average 
model ((TO]) is 

2058(s + 121200)(s + 612100)(s + 5894) 



(s 2 + 46410s + 745900000)(s + 270100)(s + 1116) 



(29) 



which has similar corresponding pole locations. The control-to-output frequency responses of the 
lifted and the average models are shown in Fig. [TBI both matching with the experimental data 
shown in Fig. 6 of |12j . 

Although the average model has matched frequency response, it still fails to predict PDB if 
v s is changed as shown next, whereas the sampled-data and the lifted models can predict PDB. 
The S plot in Fig. [T7] shows that PDB occurs if D < 0.09 (corresponding to v s = v Q /D = 25, also 
predicted by (|15p based on the harmonic balance model). 
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Figure 14: Control-to-output frequency responses; solid line for the lifted model and dashed line 
for the average model, v s = 30.84. 




Figure 15: The intersection of S(— 1, 0.4, ui p ) (solid line) and h(d) (dashed line) shows the unstable 
window of uj p £ (0.2, l)u s . 



Next, let v s = 25. The sampled-data model © is 

0.28065(z - 0.9422)(z - 0.4226)(z + 0.5305) 



(z + 1.023)(z + 0.04694)(z - 0.8816)(z - 0.9856) 

which shows PDB with two negative real poles. Lifting these two poles, the lifted model 

30795(s + 563100)(s + 101100)(s + 5961)(s 2 + 125600s + 1.318 x 10 11 ) 
(s 2 - 4558s + 9.87 x 10 10 )(s 2 + 611800s + 1.923 x 10 n )(s + 12600)(s + 1448) 

which also has a pair of unstable poles. The average model (|10p is 

17147(s + 121200)(s + 612100)(s + 5894) 
(s 2 + 303400s + 1.02 x 10 n )(s + 12810) (s + 1436) 



(30) 



is 



(31) 



(32) 
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Figure 16: Control-to-output frequency responses; solid line for the lifted model and dashed line 
for the average model. 




0.4 0.6 
Duty cycle 



Figure 17: The intersection of S(—l,D,u; s /2) (solid line) and h(d) (dashed line) shows instability 
for D < 0.09; dash-dotted line for (fl~3]) and dotted line for (fl~4"|) are approximate S plots. 

with poles in LHP and fails to predict the instability. 

With v s = 25, the S plot as a function of u p in Fig. [18] shows an unstable window for 0.36 < 
ojp/oj s < 0.54. To avoid the unstable window, the S plot shows the required slope is h(d) = 185000 
(or V h = 1.85). 

Next, let h(d) = 185000, the sampled-data poles are -0.982, -0.049, 0.881, and 0.986, indicating 
that PDB is avoided as expected. The S plot also shows that if u> p /uj s > 0.54 or uj p /lj s < 0.36, PDB 
is avoided by operating outside the unstable window. Let oj p = 0.55w s , for example, the sampled- 
data poles are -0.991, -0.036, 0.882, and 0.986, indicating that PDB is avoided as expected. □ 



6 Conclusion 
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Figure 18: The intersection of S(— 1, 0.75, uj p ) (solid line) and h(d) (dashed line) shows the unstable 
window of u p £ (0.36, 0.54)cj s . 

Sampled-data modeling and harmonic balance modeling are applied to analyze ACMC buck con- 
verters. In the sampled-data model, the system dynamics is derived directly and exactly from the 
exact switching model. The orbital nature of the nominal periodic solution is preserved. Orbital 
stability is studied and is unrelated to the ripple size of the current-loop compensator output. A 
new continuous-time model "lifted" from the sampled-data model is also derived, and has frequency 
response matched with experimental data reported previously. An unstable window of u p (current- 
loop compensator pole) which leads to PDB is found by simulations. PDB can be predicted by 
the sampled-data model, the lifted model and the harmonic balance model, but not the average 
model if the system dimension is not increased. Similar analysis can be applied to other types of 
converters. 

Accurate PDB boundary conditions, such as ([7]), ([8]), (fT2j) - (fT5j) are derived. The boundary 
conditions, which also lead to a new "S plot", greatly assist converter design to avoid instability, 
just like the popular Bode plot. The S plot has many expressions: ([7]) is for any DC-DC converters 
and is exact, ([8]) and (fT2|) - (fT5|) are for buck converters, ([8]) is exact in terms of matrices, (fT2|) is 
exact in terms of signal harmonics, and f|13|) - ()15j ) are approximate. 

With some converter parameters, the average model may show frequency responses matched 
with experimental data. The good match in one condition does not necessarily mean that the good 
match occurs in all other conditions. If one converter parameter is adjusted such that the converter 
has a negative real sampled-data pole, the average model without increased system dimension does 
not accurately predict the converter dynamics. A negative real sampled-data pole implies that the 
system is oscillatory. The traditional average model without increased system dimension does not 
have a corresponding single pole mapped from the negative real sampled-data pole. This deficiency 
is removed in the lifted model. In the lifted model, the system dimension is increased to generate 
additional degree of freedom to have better and matched frequency response. Therefore, if the 
converter does not have a negative real sampled-data pole, the average model is adequate. If the 
converter has a negative real sampled-data pole, the sampled-data and the lifted models will give 
more accurate results. 

The following design procedure for ACMC is suggested. First, use (|14p or (|15p for an initial 
design. Compared with the exact condition (|12p . the conditions (|14p and (|15p are conservative, 
but they are not as conservative as the traditional guideline (|16|). Next, use the S plot ([SD or (|12|) . 
as a function of D or oj p , to confirm that PDB is avoided. Increase the ramp slope above the S 
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plot to have sufficient stability margin. The S plot can be also used for pole assignment. If an 
unstable window of oj p exists, generally between 0.3ui s and Q.5u s , operate outside the window or 
adjust converter parameters according to ([HP , such as decrease K c R s v s , or increase oj z L. Then, 
plot the frequency response of the lifted model © to assist the voltage-loop compensation design. 
Knowledge about the sampled- data pole locations are also useful for digital control of the converter. 
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